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Abstract 



The stability of fermionic charge carriers on cosmic strings is considered. 
We show that neutral fermion currents in cosmic strings are always chiral 
or time-like, in contrast to the case of bosonic currents. The spectrum of 
bound states on an abelian SO(10) string is determined both before and after 
the electroweak phase transition. We determine the mass acquired by the 
zero mode at this transition. A range of charge carrier scattering processes 
are considered and corresponding decay rates calculated. Couplings between 
the carriers and the electroweak sector generate scattering from the plasma 
which can thermalise some currents. If the zero mode is isolated from the 
electroweak sector, then it survives. Current-current scattering is considered, 
but found to be unimportant in realistic settings where the string density is 
low. 
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I. INTRODUCTION 



Although much of the evolution of the Universe is well understood, there are still many 
cosmological phenomena for which a completely satisfactory explanation has yet to be found. 
Topological defects, such as cosmic strings, could provide mechanisms for structure forma- 
tion, CMB anisotropy, and high energy cosmic rays |1[|. Such defects form in many realistic 
unified theories. 

In the past it has been difficult to evaluate the usefulness of such ideas due to a lack 
of data. This is now changing, and predictions of CMB anisotropies from simple cosmic 
string models have been made 0. While these predictions show poor agreement with the 
observations, they do not take into account the full physics of strings models. Indeed, recent 
analysis which includes the effect of particle production as an energy loss mechanism from 
the string network show much improved agreement with data ||. One possibility is that 
the strings carry conserved currents ||. These currents will alter the evolution of a string 
network, which could lead to better agreement with observation. Indeed, an analytic analysis 
showed that a much denser string network results for electromagnetically coupled strings || . 
However, one significant implication of conserved currents is that they can stabilise loops of 
string. If persistent, these stabilised loops or 'vortons' can easily dominate the energy 
density of the Universe, placing stringent constraints on the parameters of the model. An 
analysis has been made of the implications of this for particle physics models predicting 
current carrying strings [^]. 

Fermions are a natural choice for the charge carriers of such currents. Fermion zero modes 
exist in a wide class of cosmic string models. The fermions can be excited and move along the 
string, thus resulting in a current. Consequently, fermion conductivity occurs naturally in 
many supersymmetric and grand unified theories, such as SO(10). Most attention has been 
given to massless, chiral currents, as these are most likely to be stable. In this paper space- 
and time-like currents are also considered. These naturally occur in the bosonic models of 
superconducting strings 0. In section [TI| we show that space-like fermion bound states do 
not exist in the string core in any cosmic string model. We contrast this with the situation 
for bose current carriers. In section [TIT] the spectrum of time-like currents is investigated for 
an abelian string model. We determine the complete spectrum of fermion modes, including 
both zero modes and massive bound states in our analysis. Whilst traditionally, zero modes 
are considered for current carriers, low-lying bound states can also carry currents on strings. 
This analysis is applied to the SO(10) model. The effect of the electroweak phase transition 
on fermion currents is investigated in section [TV| . 

One criticism of fermion currents is that, unlike scalar boson currents, they are not 
topologically stable. It is possible that they could decay, either directly into particles off 
the string, or through interactions with the surrounding plasma. If the decay rate is too 
high, currents will not last long enough to have any significant effect. On the other hand 
if the decay rate is too low the Universe could become vorton dominated. Such decays 
are examined in section [V]. We first examine the effect of plasma scattering on the string 
current. We show that this process can remove current carriers close to the phase transition, 
but not otherwise. We also consider the effect of current-current scattering and apply our 
analysis to a network of strings shortly after the phase transition. We show that this also 



has a negligible effect on the current stability. The results are summarised in section VI 
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II. POSITIVITY CONDITION 



In this section we consider the possible fermionic currents that can be carried by a string. 
Consider the following fermionic Lagrangian for a set of fermions Vi : 

C = ^ipA + ^mg fc ^ fc + (h. c.)] + [fc&mg^ + (h- c.)] ■ (2.1) 

— T 

Here, Vf = Ci\) i , and C is the charge conjugation matrix. 

For shorthand we write, Mf k = 0jmg fe and = 4>jmf- k . The Dirac equation is then, 

7°<^ + \l Q Mg + M fc ? 7 U + [ 7 °M^ - rC T M^C-'W k = . (2.2) 

For concreteness we work with the Dirac representation of the gamma matrices. In this 
basis 7°, 7 3 and C = ry 2 7° are real. Combining the Dirac equation and its charge conjugate 
into a single matrix equation we find, 

Ho + l°l 3 P3 + H m }(t)=0, (2.3) 



where 



4>t 



iD \ h fiD 3 



iD*J ' d V iD* 



~ (i 1 °T a D a + 1 Mg + Mj? 1 ° 1 M t % + 1 &M% T C \ 

^ m " I M^7° + &Mg*C>f i 7 + j>&M§?C + &MgC>f ) ■ {IA) 

The index a runs over the values 1 and 2. 

The operators H and P3 are clearly hermitian. If we restrict the mass matrices so that 
the Lagrangian contains only Lorentz scalar or pseudoscalar terms then H m is also hermitian, 
and{ 7 ° 7 3 ,^ m } = 0. 

Let a state |V) now represent the vector (Vi, V 4 C ) T . For a simultaneous eigenstate of H 
and P3 we have, 

(V| H H |V) = (VI ( 7 ° 7 3 P3 + P m )( 7 ° 7 3 A + #J IV) , (2.5) 

or 

w 2 = fc 2 + ^| {7 o 7 3p 3) Hm} + ^| H m H m |V) • (2.6) 

Now 

{ 7 °7 3 A, H m } = [P 3 , 7 ° 7 3 # m ] + { 7 °7 3 , • (2.7) 

For the mass matrices we are considering, the second term vanishes. Since | V) is an eigenstate 
of P 3 , 

(V| [P 3 ,7°7 3 ^ m ] \1>) = M [h^WH m ] |V) = . (2.8) 
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Thus the anticommutator term vanishes and we have 

w 2 = k 2 3 + (^\H m H m \^) . (2.9) 

As we have seen, H m is hermitian, thus the expectation value is a weighted sum of the 
squares of real eigenvalues, i.e. it is positive definite: 

w 2 > k\ . (2.10) 

Thus for scalar and pseudoscalar mass terms, we see that there are no space-like fermions 
on the string and the total energy-momentum of the charge carriers is null or time-like. We 
can interpret this result physically as guaranteeing the stability of the string against the 
spontaneous formation of a fermionic condensate. This contrasts sharply with the bosonic 
case || , where the bare string is unstable to the spontaneous formation of a bosonic conden- 
sate. This not only allows the formation of a uniform, static condensate, but also space-like 
excitations, i.e. static condensates whose phases wind along the string. The presence of such 
states relies on two features of the bosonic model which are not present in the fermionic 
case; the non-linearity of the equations of motion and the effective scalar mass squared 
being negative in the string core. 

Of course if there are oppositely charged fermion carriers moving in different directions 
along the string, the charge current can be space-like, but the energy-momentum is still null 
or time-like. 



III. BOUND STATES IN THE ABELIAN STRING MODEL 

Having shown that only light- and time-like fermions exist on cosmic strings, we will now 
determine the complete spectrum of currents in a simple U(l) string model. A model with 
a suitable symmetry breaking has the Lagrangian 

C = (D^s)W^ s - -F^ - ~ (|$ s | 2 - v 2 ) 2 (3.1) 

where -D M $s — [9^ — ieA^)Q s and F^ v = 9^A V — d u A fl . This is just the abelian Higgs model. 
In the usual vacuum solution the mass of <3>s is M s = VXr], and the gauge field has mass 
My = y/2er]. The vacuum is topologically non-trivial, allowing cosmic string solutions to 
exist. To simplify the analysis, we consider strings with winding number 1. In this case 
(Afj) = 5^a(p)/er and ($s) = vf(p) eld > where p = Mgr. The boundary conditions of the 
two functions are /(0) = a(0) = and /(oo) = o(oo) = 1. The radius of the string is of 
order Mg" 1 . 

The model can be extended by adding a Weyl fermion field, The extra fermionic part 
of the Lagrangian is then 

^fermions = ttW^tt - h 9 ^ S ^ C + Ug*^^ , (3.2) 

where = (I, a 1 ), # c = ia 2 ^* and = (<9 M - \ieA^)m . 
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The variation of <3>s means that the cosmic string acts as a potential well for the fermions. 
As well as the usual continuum states, there will be additional fermion states which exist 
only on the string. These are the fermion zero modes and massive bound states. The fermion 
field \l/ can be expressed in terms of momentum eigenstates. The bound states have only z- 
and angular momentum (k z , n) since they are restricted to the string core. Apart from the 
massive bound states there are also massless chiral solutions. These have n = and, as a 
consequence of their chirality, only travel in one direction along the string. 

* = E (cikUite- 4 ** + clUi*e iWt ) 

k z ni 

+ (cok^oke"™' + 4k^ok e ™') + E( continuum states ) • ( 3 - 3 ) 

k z >0 k 

The index % runs over the different masses of the bound states. It can be seen that the 
massless chiral states (C/ok) are rea l- The field equations for the spinors U and U can be 
found by varying ( |3.2|) with respect to The wavefunctions are then found to have the 
form, 



jj = M I Vw + k zX i{p) \ e i(k z z+ne) 
lk V^L V Ww - k zX 2(p)e ie J 



JT - M ( VwTh_X3{p) \ -i(k z z+n8) ( n A) 

^"^(J) 6 " 7 ^^' (3 - 5) 

where L is the length of the string, M _1 is the effective radius of the wavefunctions and 
M = \gf\r] ~ Ms is the fermion mass away from the string. The spinors are normalised so 
that f d 3 x (\U ik \ 2 + \Uikl 2 ) = 1 and Jd 3 x\U 0k \ 2 = 1/2. 

The allowed masses of the bound states (Mb) can be determined by finding normalisable 
solutions for the equations satisfied by the functions Xi- m a massive bound state's rest 
frame these are 

(Tl d \ 

d p-- + Yp) Xl + MbX2 + Mfxs = (3 ' 6) 

M s (d p + -)x2-M BXl -Mf X 4 = (3.7) 



P 2p, 
n a 

~P + Z~P / 
n — 1 

7~ 2p, 



M s \ d P + - + — ) X3 + M B X4 + Mfxi = (3.8) 
M s ( d p - — - ^ ] X 4 - M B X3 - Mf X 2 = Q (3.9) 



Of course n must be an integer for the solution to be single valued and hence physical. 
Examining the above equations we see that given one solution, another can be obtained by 
putting n — > —n, Xi ^ X3 an d X2 ^ Xi- Thus the solutions will occur in pairs. 
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FIG. 1. Spectrum of n = 0, ±1 massive fermion bound states. 



We will use a variation of the shooting method to determine the allowed values of Mb. 
At large p the solutions of ( |3.6|) -( |3~9D have exponential behaviour. Two of them decay and 
so are acceptable. In the case of the small p solutions, only two of them give a normalisable 
state. 

Each of these 4 solutions can be numerically extended to some intermediate value of p 
(of order 1). It can then be seen if any non-trivial combinations of the large and small p 
solutions match up there. We find that the values of Mb for which normalisable solutions 
exist satisfy 

M 2 - M 2 = (M + a ni M s ) 2 > (3.10) 

where the a n i are functions of M/Mg. Figure [I] shows the variation of a n i and the number of 
bound states with respect to the fermion to Higgs mass ratio (M/Ms). Each line corresponds 
to two bound states. For simplicity we have taken the Higgs and gauge field masses (Mg 
and My) to be equal. With this choice of parameters all solutions have either \i or X4 
identically zero. Equations (|3.6| )- (|3.9| ) can then be reduced to a second order problem. The 
solutions with % 4 = are shown in the figures. It can be seen from figure p] that the number 
of bound states increases with the size of the off-string fermion mass. For small values of M 
there are just two massive bound states in addition to the chiral zero mode. 

We find that in the region of parameter space in which a ni is roughly constant, M 2 m 
M$JM 2 — Mf>/2. For the massless states, when Mg = My, this expression for M is exact 
and can be proved analytically. In this case the two sides of ( |3.5| ) are actually equal. 
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2. Zero mode and lowest massive bound state solutions in the abelian string model with 
: 2. 



Plots of the solutions for M/Ms = 2 are shown in figure [| As expected they decay 
outside the string. We have thus found the full spectrum of fermion bound states for this 
model. 



IV. CURRENTS AFTER THE ELECTROWEAK PHASE TRANSITION 

The toy model of the previous section can be embedded in a phenomologically credible 
grand unified theory, such as SO(10). Two suitable symmetry breakings which can give rise 
to cosmic strings are 

SO(10) SU(5) x Z 2 — > Q SM xZ 2 ^ SU(3) C x U(1) Q x Z 2 , (4.1) 

SO(10) - SU(5) x U(l) - Gsm x U(l) ^ Gsm xZ 2 ^U SU(3) C x U(1) Q x Z 2 (4.2) 

with Gsm = SU(3) C x SU(2)l x U(l)y. The discrete Z 2 symmetry allows the formation 
of topologically stable cosmic strings. In this case the string gauge field is a neutral GUT 
boson, and $s transforms under the 126 representation of SO(10). The electroweak Higgs 
field $ transforms under the 10 representation. 

The fermion sector of the SO(10) GUT contains all the Standard Model fermions, and 
an extra right-handed neutrino (z/ c ) for each family. Only right-handed neutrinos couple to 
($s)j while neutrinos of either helicity couple to (<&). 

The abelian string's gauge field has a non-trivial effect on the electroweak Higgs field. 
The components of $ have charges ±1/5 with respect to the generator of the GUT string 
gauge field, so ($) will not be constant in the presence of a string. A non-zero Z field is also 
required to give a finite energy solution fllPf . 

The r-neutrino mass terms in a string background are 

(^(xt) m "o w )(:r w> 
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m ew ~ 10 2 GeV is the electroweak energy scale and M ~ 10 16 GeV is the GUT energy scale. 
For simplicity we will neglect mixing between different families. The function h gives the 
radial dependence of the component of ($) which couples to the neutrinos. The electroweak 
gauge field has (Zg) = — 6(p)/(v40er). The boundary conditions of h and b are 6(0) = and 
h(oo) = b(oo) = 1. Inside the string h « (m ew /M)c =constant, with c ~ 10 4 |p.l| . While 
the GUT fields take their vacuum expectation values outside the string (whose radius is of 
order M _1 ), the electroweak fields vary over a far larger region with radius of order m~^. 

Since e = m ew /M <C 1 the neutrino mass eigenvalues outside the string (or in a vacuum) 
are approximately M and = m 2 ev /M. The corresponding mass eigenstates are then 
approximately v c + ev and v — ev c . This illustrates the seesaw mechanism |[L2|| . Although 
the r-neutrinos have the same couplings to the electroweak Higgs field as the top quark, 
the GUT Higgs ensures that v c is superheavy and v is very light, as is required to agree 
with observation. Recent measurements have suggested that v does indeed have a small 
mass 
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Once $ gains a non-zero expectation value, the v c currents are no longer solutions of the 
field equations. We will start by considering the fate of the massless v c currents. It seems 
likely that they mix with a v state to give a low mass bound state. We denote the resulting 
mass of the bound state by pM, where the dimensionless parameter /i is expected to be very 
small. Using the ansatz 

; ,c _ ( Xi(p) cos(pMt) \ ( Ci(p)e- ie cos(pMt) \ 

[ X2 (p)e ie sm(pMt) ) ' ^ Up) sin(pMt) J 

we can look for solutions which reduce to the zero mode solutions of ( |3.6| )- (|3T9| ) when ($) = 0. 
Putting Ms = M for simplicity, the field equations reduce to 




Xi + PX2 + ehCi 


= 


(4.5) 


X2 - PXi - 


= 


(4.6) 


Ci + K2 + ehxi 


= 


(4.7) 


C2 - pCi ~ thx2 


= Q 


(4.8) 



As with the abelian string, the existence of bound states can be investigated by examining 
the large and small p solutions of ( (4.5| )-( ^8| ), and then trying to match them at intermediate 
p. There are two well behaved solutions at small p and two at large p if /1 < e 2 . 

We will attempt to do this by finding approximate solutions in a simple 'top hat' ap- 
proximation of the string background. This has f(p) = a(p) = 6(p), b(p) = 6(ep) and 
h(p) = ec+ (1 — ec)O(ep), where 9 is the Heaviside step function. It is then possible to find 
analytic solutions for p > e _1 . For smaller p the solutions can be expanded in terms of the 
tiny parameter e. 

Requiring the approximate solutions to match up at p = 1 and p = e _1 will give an 
expression for p. It is satisfied by p, — \f2^ce l& ^ with 7 = er(6/5, 1) — 1/2 ps 1.6. V is the 
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incomplete gamma function. To leading order in e, the corresponding approximate solution 
is 



X 2 = ^(7-I)e 16/5 P 



C 2 

Ci = - r P 

( 2 = -V2e 6/5 (4.9) 



for p < 1, 




~7F. 



for 1 < p < e 1 , and 




(4.10) 



^V^, f^<^o ( r'/^ (4.11) 



for p > e -1 . is a confluent hypergeometric function. 

The mass of the new bound state is Mp ~ 10~ 16 eV. Unlike the original zero mode, it is 
a mixture of left- and right-handed neutrinos. Similarly, the massless v c and z/i currents will 
also gain tiny masses at the electroweak phase transition. Since these states now couple to 
light off-string states, they can escape from curved strings ||. Applying the arguments in 
ref. || the rate of decay of states with energy of order Ms is calculated to give a lifetime of 
order 10 _40 s. This would prevent the Universe becoming vorton dominated. In addition the 
electroweak bound states are spread over a far greater region than the original v c current, so 
interactions with other fields will also be increased. If any massive bound states survive the 
phase transition they will also rapidly decay into left-handed neutrinos, although it seems 
more likely that they will simply be changed into free states by the transition. 

The above mechanism will work in other theories where zero modes mix with free massless 
particles. If there is mixing between two massive species of right and left moving currents, 
zero modes can mix to give a low mass bound state. However, as there is no low mass 
off-string state these cannot escape. Such theories always feature zero modes moving in 
opposite directions. These can interact with each other, which reduces the current. It is 
still possible for a net current to persist, and so fermion currents may provide constraints 
on these models after all. 

If a zero mode is isolated from the electroweak sector, then it persists. Consequently 
the resulting currents also persist. In particular, zero modes in grand unified theories or 
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" (b) " - (d) 

FIG. 3. Feynman diagrams contributing to plasma scattering (a,b), massive bound state decay 
(c) and vorton decay (d). 



super symmetric theories where there is no coupling between the zero mode and the elec- 
troweak sector, result in currents which are not destroyed. Generally, such currents give rise 
to stable vortons, which constrain the underlying particle physics theory ||. This is likely 
to be particularly true in D-type supersymmetric theories. 



Having established the spectrum of fermion states in an abelian string background, we 
now turn our attention to the stability of charge carriers on these strings. In the absence 
of other particles, currents carried by zero modes on isolated, straight strings are stable on 
grounds of energy and momentum conservation. However, in a realistic setting there are 
many processes which can depopulate zero modes. Strings are not isolated and in the early 
Universe they and their bound states will interact with the hot plasma. Also, bound states 
on different strings, or different parts of a single curved string, can scatter from one another. 
If there are couplings within the theory that allow a heavy neutrino zero mode to scatter 
from a light plasma particle to produce a light fermion-antifermion pair via an intermediate 
electroweak Higgs boson, we can have charge carrier decay from the Feynman diagrams in 
figures 0a and 0b. Also of interest is the decay of massive bound modes into light fermions 
(see figure 0c). In theories which allow such interactions, the charge carriers can decay. It is 
also possible for heavy neutrino currents on different strings to decay by exchanging a Higgs 
particle, as in figure |3]d. The aim of this section is to calculate these decay rates. 

The details of the calculations are given in the appendices, here we highlight the impor- 
tant differences between the calculation in the background of the string and the correspond- 
ing trivial background case. We initially consider our system to be restricted to some box of 
finite volume, V. Following the canonical procedure we decompose the field operators into a 
sum over orthonormal wavefunctions and creation/annihilation operators. Throughout we 
normalise our wavefunctions according to, 



V. DECAY RATES OF CHARGE CARRIERS 





We take the following decomposition for scalar fields, 



i 



(a k k (x)e- iurf + 6 k 0_ k (x)e™*) 



(5.2) 



k 



V2w 
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where w 2 = k-k + m| fc . The necessary commutation relations are, [ok, a k /] = [b^, b\,,] = 5^. 
The corresponding decomposition of fermionic fields is, 

*(t, x) = ]T (ckC/k(x)e-^ + 4\/ k (x)e rat ) , (5.3) 

k 

where the states and Vk are spinor valued and the sum is over momentum states. In this 
case we impose anticommutation relations {ck, c k /} = {d k , d\,} = 5 k k'- 

With these normalisations we have a simple interpretation of the amplitude, 

A = (1,21513, 4) . (5.4) 

The probability of the interaction, characterised by S, converting the initial state, |1,2), 
into the final state, |3,4), is simply, 

P=|-4| 2 - (5.5) 

We can now consider a simple, second order tree diagram in $\I/t\|/ c theory, where the 
interaction is given by, 

H mt = ig<S>^\ia 2 y* 2 + ig*^*^ia 2 ^ 4 + (h. c.) . (5.6) 

The amplitude is given by, 

A = (in| T\g\ 2 i 2 J d 4 x^\(x)ia 2 ^ 2 (x)^j(x) J d 4 y^(y)ia 2 ^ 4 (y)^(y) M> , (5.7) 
where incoming and outgoing states have the form, 

(in| = (0|cica , |out) =c!4|0) . (5.8) 
Expanding the field operators in the usual way we obtain, 

A = M'E / d " x J d 4 x' Ul(K)ia 2 U;(K)e i ^ +w ^ t 

k 

x (0| T^x^y) |0> Uj(x')ia 2 U 4 (x')e- i(ws+W4)t ' . (5.9) 
Expressing the Green's function as a sum over a complete set of states, we have, 



A=\g\ 2 Y, J d 4 x J d\r'rUx)icr 2 r.Ux)<' 

k 



x tt i;(X 2 2 1 WU^)e~^^' , (5.10) 

where T is the temporal extent of the region we are considering and is the (constant) 
effective mass of the scalar mode. 

It is useful to consider briefly how this calculation would proceed in a trivial background. 
In a space-time box of volume VT, we have, 
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k 



(2tt)< 



Each wavefunction introduces a factor of V x l 2 to the amplitude. The spatial integrations 
yield finite volume approximations to energy and momentum conserving delta functions 
which we denote b\ v , where J e ik x d 4 x = (2ir) 4 6$- v (k) and 5^(0) = TV{2ix)^. 

The value of the propagator is fixed by energy /momentum conservation and, up to 
dimensionless factors, we have, 

A ~ l^^-S^ih +k 2 -k 3 - k 4 ) , (5.12) 

where ki = k\ + A; 2 = h$ + £4. 

The interaction cross section is given by, 



final states 



Replacing this sum by integrals over the final state momenta and including all the spinor 
factors, we find, 

\n\ A k 2 

- - 151 kl (5.14) 



16-7T (kj — rri^) 2 

Now we can consider a situation of interest. If the incoming particles are a zero mass 
particle of energy w zm and a light plasma particle of energy of order the temperature, then 
kj ~ w zm T and, 

\g\ 4 w zm T 

& ~ ~, ^ ^T77 • (5.15) 

(w zm T - m\) 2 y J 

Having discussed the calculation in a trivial background, we can now consider the string 
background. The major difference arises in the spatial integration associated with the posi- 
tion of the initial vertex. Let particle 1 be the bound state, 

~ jmt-ikuz-imo e -Mr ^ (5 16 ) 

Vl 

where M~ l is the radius of the bound state. As this state is localised close to the string, 
it will only overlap significantly with the lowest angular mode components of the incoming 
scattering wave and so the integration over the position of the initial vertex will be dominated 
by these modes (and for the same reason the lowest modes of the propagator). 

We consider the incoming particle to be a plane wave asymptotically and do a mode by 
mode matching of this incoming plane wave onto the states in the string background. As 
we only have square integrable wavefunctions in the string background, in general we must 
also include an outgoing scattered wave: this is the source of Aharanov-Bohm scattering. 
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The general features of both the incoming particle and intermediate particle wavefunc- 
tions are a growth region for k^r ^ 1 and an oscillatory region for /c^r •> 1, where k^ is the 
momentum of the particle transverse to the string pi. The transverse momentum of the 



incoming particle will be much less than M, so we expect some power law type behaviour 
in the region of overlap with the bound state. Similarly, if the transverse momentum of 
the intermediate particle, kjT is smaller than M, its wavefunction will also be roughly a 
power law in the important overlap region. The integration over the initial vertex position 
then produces a number that is suppressed due to the limited extent of the bound state 
wavefunction, but is only mildly dependent on kj?- Only when kiT exceeds M do we find 
a wavefunction oscillating in the overlap region. Thus only for kjT ^ M does the integral 
decrease. Instead of a momentum conserving delta function, this spatial integration yields 
something more like a step function, which permits transverse momentum non-conservation 
up to the mass scale M. This behaviour is natural as the string breaks transverse translation 
invariance and the object bound to the string has an effective radius M~ l . 

The violation of momentum conservation in the string background opens up significant 
areas of phase space that are forbidden in the trivial background. Of particular importance 
is the possibility of resonant scattering. We consider first the decay of a charge carrier 
into a fermion and a Higgs boson. In the string background transverse momentum can 
be acquired from the string to place the Higgs particle on shell in a much larger region of 
phase space. This introduces a factor of M 2 M B . The small size of the region in which the 
wavefunctions overlap gives an extra factor of M~ 2 . The full decay into three fermions can 
then be considered to be an initial decay into a light fermion and physical electroweak Higgs 
boson, followed by Higgs decay. In appendix |CJ it is shown that the massive bound state 
lifetime is 

r~(\gA 2 M B y l , (5.17) 

where g u is the Yukawa coupling in the neutrino electroweak mass term. Thus (|5.17|) results 
in a small lifetime. 

In arriving at this lifetime we have made certain assumptions, in particular we have 
neglected the possibility that the incoming and outgoing fermionic wavefunctions may be 



amplified near the string ||14|| . We have also neglected back- reaction on the string. As 
a significant contribution to the amplitude comes from large transverse momentum non- 
conservation, the use of a static string background might not be a good approximation. 
However, if the string recoils it will absorb some of the energy of the interaction. This will 
reduce the amount of momentum non-conservation required to make the outgoing particles 
on-shell. Thus these approximations are not likely to increase the lifetime, and ( |5.17| ) should 
therefore be taken as an overestimate. 

The fate of massless bound states is also complicated by momentum non-conservation. 
Of particular interest are the bound states that stabilise vortons. Contraction of the vorton 
will ensure that the states at the Fermi surface will have GUT scale momenta. The lifetime 
of these high momentum states is critical, if they decay the vorton will contract and promote 
low momentum states to high momentum. Energy and z-momentum conservation prevents 
the massless states decaying spontaneously, in contrast to the massive bound states. It 
is however possible for them to decay by interaction with plasma particles or other zero 
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modes. If Ms > M it is also possible for high energy massless currents on a curved string to 
decay by tunneling to free heavy neutrinos. However the rate will not be significant unless 
M S »M| 

If we consider one of these states scattering from a typical plasma particle, we have 
a centre of mass energy of order \/M$T, well above the mass of the intermediate parti- 
cle. Transverse momentum non-conservation again allows for resonant scattering. Including 
amplification of the incoming plasma particle wavefunction, the lifetime of these high mo- 
mentum zero modes is found to be, 

M 2 (T\ 2Q , 

) , (5-18) 



where Q is the charge of the plasma particle under the string gauge field. We have taken the 
plasma particles to be massless and ignored any temperature dependent corrections to their 
mass. Providing ^w zm T > ~ T, resonant scattering is possible and the cross-section is 
largely independent of the zero mode energy, w zm . The unamplified lifetime of the modes 
thus scales with the plasma density. Conversely, the amplification factor decreases with 
increasing temperature as the ratio of the GUT scale to the typical thermal energy grows. 
In the radiation dominated era we can take t = aT~ 2 and the lifetime becomes, 

r ~ M 2 - 2Q \g u \- 2 T 2Q - 3 = M 2 - 2Q \g v \- 2 a Q - 3 l 2 t 3 / 2 - Q . (5.19) 

The probability of a zero mode state scattering in time interval dt is dt/r, thus the probability 
of a zero mode state scattering after some time ti is, 

-n*± f ~ , t Q - 1/2 1 

P(decay after U) = 1 - e Jt > T = 1 - exp I -M 2Q - 2 \g v \ 2 a 3/2 ' Q ^ — - \ . (5.20) 

If the magnitude of the exponent is small there is a small probability, thus zero modes are 
stable if, 

M 2Q ~ 2 1 g v | 2 a 3/2 ~ Q l 2/(1 ~ 2<3) \ . (5.21) 



U > O 

Now, a ~ Mpi/10 and in the case of SO(10), Q = 3/10, leading to the condition, 



U > O ([M-%| 10 a 6 ]) > O 



M P r" 



10M 



9A 10 M- X \ . (5.22) 



As the lifetime varies only slightly faster than T~ 2 , this result for £, is very sensitive 
to the Yukawa coupling. For M ~ 10 15 GeV, if \g u \ = 1 zero mode states populated after 
ti ~ 10 15 t GUT will be stable, while if \g u \ < 0.03, this scattering is never significant. In the 
SO(10) model g v is also the Yukawa coupling for the corresponding quarks, thus there is an 
epoch when z/: zero modes will scatter from the string, but v c e and v c zero modes will never 
scatter by this process. Thus the interaction with plasma particles can not significantly 
remove zero modes from the string. Note that it is also possible to create currents using the 
above interactions in reverse. Hence, if thermal equilibrium is reached the number density 
of zero modes will be of order T. 
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Within the SO(10) model there is also the possibility of mediating these processes by 
GUT mass Higgs fields with zero VEV. In this case the Yukawa coupling need not be small, 
but the centre of mass energy of the interaction is only of order the intermediate particle 
mass for T ~ T GUT . Thus below the GUT temperature the reaction rates for these processes 
are rapidly suppressed by powers of T/T GUT . 

None of these plasma scattering processes can remove v c e and v c zero modes, and are only 
significant for v° immediately after the phase transition. Thus, they are unable to prevent 
the vorton density from dominating the energy density of the Universe. 

The plasma scattering processes considered above failed to remove zero modes due to 
the decreasing plasma density at late times. A distinct category of process is the scattering 
of a zero mode on one string by a second zero mode on another string. This is particularly 
relevant for vortons as they form small loops with a typical radius only one or two orders 
of magnitude larger than the string width. We thus have zero modes moving in opposite 
directions on opposite sides of the vorton. For simplicity the decay rate can be calculated 
by considering two straight, anti-parallel strings (appendix jD|). Physically one expects a 
suppression due to the finite separation of the initial and final vertices. As the initial and 
final vertices are confined to different pieces of string, the amplitude contains a factor of e lkTR , 
where k^ is the transverse momentum of the intermediate particle and 2R is the separation 
of the strings. When the integration over the intermediate momenta is performed, this factor 
can cause the integrand to oscillate, leading to a suppression of the amplitude. This occurs 
if kxR is large where the standard propagator factor is significant. 

If the intermediate particle is an electroweak Higgs boson, then the standard propagator 
factor peaks for kx around the electroweak scale, krR is tiny in the important region, there 
is no oscillation and no suppression. The cross-section in this case is found to be, 

" ~ w£w • (5 - 23) 

This cross-section is dimensionless as the scattering is effectively in one spatial dimension. 
In this case resonant scattering is not possible, thus the Higgs width does not enter and the 
rate therefore contains a factor of g„. 

Conversely, for a GUT scale mass as the intermediate particle, the standard propagator 
factor peaks for ~ M GUT , giving k?R ~ 10 — 100. In this regime the reaction rate is 
found to display the exponential suppression expected on physical grounds, 

The exponential suppression in Q5.24Q makes such processes irrelevant in all physical situa- 
tions. 

Taken at face value, if electroweak particles mediate current-current scattering on dif- 
ferent segments of string then ( |5.2[j| ) gives a short lifetime for charge carriers on a vor- 
ton. However, there are complications in directly applying this result derived for straight 
strings to curved vortons. The presence of charge carriers on the vorton loops results in the 
loop carrying angular momentum, which must be conserved in any physical process under 
investigation. The above calculation does not take into account conservation of angular 
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momentum, as the system considered has no rotational symmetry. However, the angular 
momentum gives rise to a centrigual barrier, suppressing the above decay channel. Com- 
bined with energy conservation, this would prevent massless modes on the string scattering 
into massive modes. However, the fermionic spectrum has not been calculated for a circular 
loop and there is no reason to expect the zero modes to remain massless and thus, angular 
momentum conservation can only be considered in a consistent rotationally invariant cal- 
culation. The calculation above works consistently with straight strings, thus while ( |5.23| ) 
may not be directly applicable to vorton decay, it is relevant for interactions on non-circular 
loops and scattering of currents on a string network. 

At very early times in the friction dominated regime, the string correlation length is 
small. The interaction of (|5.23|) could be operative. However the string density drops too 
quickly for it to be significant. As the density of current carriers builds up, the inter-string 
separation increases, so the cross-section decreases significantly. Consequently, the current- 
current interaction is unable to reduce the number of current carriers on the string and 
prevent vorton formation. 



VI. SUMMARY 

In this paper the existence and form of fermion bound states and their corresponding 
currents on cosmic strings were investigated. Only time-like and light-like fermions can 
exist on the string, in contrast to the bosonic case. Using numerical methods, the discrete 
spectrum of states for an abelian string model with one fermion field was determined. We 
found that the number of bound states increased as the Yukawa coupling grew. For very 
low values there are just two massive bound states and a zero mode. 

Bound mode states will always occur. The occupancy of these states depends on the 
decay modes of the carriers. Since the cosmic string breaks Lorentz invariance, transverse 
momentum is not conserved in interactions in a string background. This leads to an enlarged 
phase space and increased cross sections. On the other hand, the fact that the bound states 
are confined to the string reduces the overlap of the particle wavef unctions, tending to reduce 
the cross sections. 

Unless the massive carriers are isolated from the electroweak sector, they will decay and 
the states will empty on a time scale of 10 _33 s. If the carriers are stable, they will persist 
on the string and carry angular momentum, contributing to vorton formation. 

In order for the massless current carriers to decay they must interact with other particles, 
such as carriers on other strings or plasma particles. The most significant decays involve 
electroweak Higgs intermediate states. Non- conservation of momentum means that plasma 
scattering is usually resonant. However the small size of the Yukawa couplings and the 
plasma density mean that the rate of this decay is too small to be cosmologically significant 
(at least for v eR and u^ R current carriers). Decays involving GUT boson intermediates are 
also possible, although the momentum non-conservation is not large enough to make these 
decays resonant. The interaction rate in this case is also small. 

Massless current carriers can also decay by scattering with currents on other strings. The 
cross section for this interaction is tiny unless the two strings are very close together. For 
GUT mass intermediate states appropriate densities are never physically realised, but for 
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electroweak intermediate states suitable densities arise immediately after string formation. 
A situation similar to this occurs when the current is on a string loop. In this case the 
current will decay rapidly. Unfortunately this result cannot be applied directly to circular 
loops, where angular momentum conservation must be taken into account. Consequently, 
it can not be applied directly to the case of cosmic vortons. This situation requires more 
detailed analysis flSfl . 

In the SO (10) model, after the electroweak phase transition, the situation changes dra- 
matically. The right-handed neutrino zero modes mix with left-handed neutrinos allowing 
currents on curved strings to decay by tunneling into free left-handed neutrinos. This averts 
any cosmological disaster. 

However, if the zero mode is isolated from the electroweak sector, then the decay processes 
considered in this paper are not operative. As a consequence, the zero mode survives. For 
GUT scale strings, this results in a cosmological disaster as discussed in [|J. 
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APPENDIX A: BASICS: PROPAGATORS, INTERACTIONS AND 

WAVEFUNCTIONS 

1. The Scalar Propagator 

We will begin by finding the propagator, G(x, y), which satisfies, 

OG(x, y) = (d 2 - V 2 + V(r) 2 )G(x, y) = -i5\x - y), (Al) 

where 0$ = is the equation of motion. V(r) 2 represents the position dependent potential 
that the scalar field experiences in the string background. The standard connection with 
the time ordered product holds here too; 

6 (0| T<$>(x)<$>*{x') |0> = -i5\x - y). (A2) 

We now have to find a suitable representation of this Green's function. Consider the quantity, 

g — iwr 

G(x -x') = i£ 7 =7— 2 — 2 k (x)0*(x'), (A3) 
k i [w — k — m^ k \ 

where r = t — t' and the sum extends over states of all masses, i.e. we include off-shell 
states in the sum. We will work in a space-time box and so impose boundary conditions 
at t = — T/2 and t = T/2. The allowed energies are then quantised in just the same 
way as the momenta. Normalising these states to unity introduces a factor of 1/y/T into 



16 



each wavefunction. We will display this factor explicitly so that only the 1 / yV factors are 
implicit in any wavefunction. 

The full wavef unctions, <pk{ x ) — e ~ lwt( P]s.( x ), satisfy, 

6<p k (x) = \-w 2 - V 2 + V r (r) 2 ] e ~™Vk(x) 

= [- W 2 + k 2 + mJ ifc ]e-™Vk(x), (A4) 

where m^^ is the (constant) effective mass of the scalar mode of momentum k. 

To see that we have a Green's function, act with the operator on G and use the fact that 
the off shell states are a complete set of solutions of the 4-dimensional eigenvalue problem. 
We tame the singularity in the usual way with an ie to produce the Feynman propagator. 



2. The Interaction 

Now consider the simple, second order tree diagram in a theory with 

Him = ig$y\ia 2 y* 2 + ig*$*^ia 2 ty 4 + (h. c.) . (A5) 
We will start by considering incoming and outgoing states of the form 

(in| = (0| 6ic 2 , |out) = c\c\ |0) . (A6) 



Substituting in the expressions for the external fields (|5.3| ) we arrive at the VEV of the time 
ordered product in the usual way, 

A = M'E / d * x J d 4 x'ul(x)ia 2 U^)e i ^ +W2 ~ w ^ t 

k 

x ul^)ta 2 U^)e^ w ^- w ^' . (A7) 

T[kj — m^ k \ 

3. The Amplitude 

First we will evaluate the amplitude (|A7|) in a trivial background, using cartesian coor- 
dinates. We will be mainly interested in 2-component left-handed massless Weyl fermions, 
with Lagrangian £ = ^Ha^D^ . The wavefunctions in (|5.2| ) and (|5.3|) are then 



e ik-x 



'V 

Jkx ikx 



6 > 
yV yV 



with n k =-L( ) , (A8) 



where k x + ik y = k^e 1 ^ . Since we are considering massless fermions w = |k|. The sum over 
intermediate states can be approximated by an integral, 
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VT r 

E-(^/^. (A9) 

The amplitude then becomes, 

„ ( J,x-(k\+k2— ki) g-ix'-(k3+k4-ki) 

A = ^ \ 27T yv 2 J d4k i d4xd4x ' u li ia2 <2 kj-ml «k> 2 «k* • ( A1 °) 

We will define (2Tr) 4: 5^- v (k) = J e lk ' x d 4 x to be the equivalents of the usual delta functions 
at finite time and volume, with S^- v (0) = TV(2tt)~ 4: . The large time and volume limits are 
implied. 

Evaluating the integrals gives 

A = +h _ h _ h) , (A11) 

V 2 kf — 

with kj = ki + k 2 = k 3 + k 4 . 

The total cross section for the interaction is obtained by summing over the possible final 

states 

final states reI 



The sum over final states can be replaced by integrals over the final state momenta, 

V fVd'h f Vd 3 k 4 \g\\2n)\ 4 4 , J^^V/K^VJ 2 
a= T^J^rJ ( 27 r)3 V 4 ^(O)^rv(fci + - *3 - **) 

(A13) 

The relative velocity of the incoming particles is i> re i = &i /^/(wiu^). Using this and 
|w£W 2 -u k /| 2 = k fl k' fl /(2ww') = (k + fc') 2 /(4ww') gives 



/ + h - h) . (A14) 



8(2tt) 2 (A; 2 - m 2 ) 2 J w 3 Wa 

The integral is Lorentz invariant, and can easily be evaluated in the centre of momentum 
frame to give 

° = \^{k}-m\y (A15) 



4. The Higgs Width 

The part of the electroweak Higgs which couples to neutrinos can decay into two quarks, 
$ — > q + q c , where q — u,c, t. This is a first order process with rate given by 
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V 2 d 3 k l d 3 k 2 
T(2vr) 6 



z ff *(27T) 4 ^ Uk2 6(m, -wi- w 2 )5 3 (k 1 + k 2 ) 
2m^V 3 



3 E 

o V- \9 q \ 2 m^ fd 3 k 1 d 3 k 2 . , 3 , , 

3 E q/o„\2 / ~ — — o(m0-iwi-ty 2 )<) (k 1 + k 2 ) 



8(2vr)2 7 Wl ^2 



3 E kl 



q=u,c,t 



2 

16tt 



(A16) 



5. Wavefunctions in a Cosmic String Background 

Having established our normalisations and conventions in the familiar trivial background, 
we now turn to the calculation in the background of a cosmic string. We work in a cylindrical 
box of length L and volume V. The normalised on-shell states are given in section |T|. 

In the string background it is natural to perform an expansion in angular mode number, 
yet we still wish to consider an incoming plane wave. The two expansions can be matched 
by first expanding the plane wave in terms of Bessel functions; 



e 



ik(x cos j3+y sin 0) _ i(k x x+k y y) _ -q j (t. r \Jq{S-0) 



J2 i q J q (kr)e iq{0 - p) , (A17) 



and then asymptotically matching these Bessel functions mode by mode to the modes in 
the string background. As we only have square integrable wavefunctions in the string back- 
ground, in general we must also include an outgoing scattered wave. This outgoing wave 
is the source of Aharanov-Bohm scattering. Using an angular mode decomposition of the 
intermediate particle, all the wavefunctions at the initial vertex are in the form of cylindrical 
modes. 

We must now consider the participating wavefunctions. Let particle 1 be the bound 

state, 

^i(x) ~ eim.t-ik^z-inie e -Mr _ ( A1 g) 



As this state is bound to the string, it will only overlap significantly with the lowest 
angular mode components of the incoming scattering wave, and so the integration over the 
position of the initial vertex will be dominated by these components. Similarly, we need 
only consider the lowest angular modes of the propagator. 

The modes of a generic fermion wavef unction in the string background have three impor- 
tant regions. If the transverse momentum is kx, these regions are fc^r >l,l/£/r>r>l/M 
and r < 1/M. In the large r region, fc^r > 1, we have some Bessel function at large 
argument. This region dominates the normalisation integral. In the intermediate region, 
~ r ~ 1/M we have some Bessel function at small argument. The order of this Bessel 
function is shifted by the string gauge field Hill . In the lowest angular momentum modes 



there is a mixture of Bessel functions of the first and second type, leading to one component 
of the spinor varying like r~® in this region, where Q is the magnitude of the charge of 
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the fermion under the string gauge field. Finally, inside the string, r < 1/M, the spinor 
components tend to constants as r — > 0. The effect of the intermediate region is thus to 
amplify the wavefunctions by 0([M /kxy)- In the following appendices we calculate unam- 
plified cross-sections. The amplification factors, which depend on the fermion charges, will 
be added at the end. 

For the outgoing states a construction as above is less intuitive. We can think of the 
incoming wave undergoing classical scattering and some small part of it also partaking in 
the quantum scattering. The corresponding interpretation for the outgoing wave is that 
the quantum scattering excites states consisting of an outgoing plane wave and outgoing 
Aharanov-Bohm scattered waves. Whilst this does not lend itself to a clean interpretation 
in terms of two-to-two scattering, it is just another manifestation of the inappropriateness 
of the plane wave as the asymptotic state in the background of the string. The 1/r fall off 
of the gauge field of the string introduces long range interactions in the string background 
which we will neglect. This corresponds to assuming that the outgoing states are free. 

The outgoing states we are interested in are light and have little interaction with the 
string. Thus we make the approximation that the outgoing states do not interact with the 
string fields and the plane wave expansion employed for a free plane wave can be used for 
the outgoing states. 



APPENDIX B: SCATTERING FROM A MASSLESS BOUND STATE ON A 

STRING 

We now consider the problem of real physical interest: a plasma fermion scattering from 
a bound fermion via an intermediate scalar to produce two light fermions. We consider 
theories with couplings of the form: 

H mt = ig v ^Ua 2 % + ig u ^ u ia 2 % c + (h. c.) , (Bl) 

where \lV is the heavy neutrino, $ is a light Higgs field, and the other \l/ are light fermions. 
In the SO (10) model the couplings g v and g u are the same. 

For the SO(10) model the scattering of interest is u^ tT + v — > u + u c , where v^ tr is 
the massless current carrier. Figure [3]a shows the corresponding Feynman diagram. The 
scattering amplitude is then 

A = (u s c tr (k 1 )u(k 2 )\S\u(k 3 )u c (k 4 )) 

k 

xG{x - ;rO^£(xV 2 ^k 4 (xOe- 4{ ™ 3+W4)t ' . (B2) 

To keep things simple, we will use plane wave approximations of the Higgs and light fermion 
states (|A8|) . Using these and the expression for £/ok (|3.5|) , the above amplitude can be 
expanded to give a similar expression to ( |A10| ). Most of the integrals can be done in the 
same way as before, leading to 
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A = ^:- 7 =£=(2tt) 2 | d 2 x5^\h + k 2 -k 3 - h) 

X W^V / 1 _^2 e -^ e _ r M/ + ^ (B3) 

kj -m^ + im^T^ V ^2 

where k • x^ 2 -* = k x x + k y y. The Higgs particle momentum satisfies kj = k 3 + k 4 and 
{ki) z ,t — {k\ + k 2 ) Zj f The above <i 2 x integral is most easily evaluated in polar coordinates. 
Using a Bessel function expansion of a plane wave ( |A17|) gives 

f e i(k x x+k yy ) d0 = 27l J (kr) (B4) 
Using a numerical solution for the r dependence of the zero mode, we find 

J (kr)e-f M f+^ dr 'rdr « ^e- fc2 /( 2 ^ . (B5) 



Co ~ 0.7 is a slowly varying function of M/M$. These two results allow the integrals in ( p3| 
to be evaluated, giving 



A = g u g:-^=^\ l - — e-^5^\k x + k 
u V2ttLV 3 MV w 2 



2 



x 9 <^f u ^ c -(fc,-fa)|/<M?») . (B6 ) 

kf-m^ + im^T^ 

The corresponding expression for scattering in a trivial background has a delta function 
instead of the gaussian. The physical significance of this is that transverse momentum is 
not conserved at the first vertex. This is to be expected since the string breaks transverse 
Lorentz invariance. The extra momentum is absorbed or provided by the string, up to the 
scale M = y/M s M /2. This is the effective energy scale of the string with respect to bound 
state interactions. There is also an extra M~ l factor, which is a combination of the bound 
state normalisation and the fact the interaction at the first vertex is confined to the string 
core, whose effective area is M~ 2 . 

The total cross section is again obtained by squaring the modulus of the amplitude and 
integrating over the final state momenta. After squaring it is convenient to introduce two 
extra ^-functions 



UP _ | fl a 12 T ( 27 3c _ **A k J c -fe-fe)^/M 2 

l9u9ul V* AM 2 { w 2 ) Ikj-ml + tm^^ 6 

x f f^Ls\k! -k 3 - h) . (B7) 
J W3W4 

In this case t> re i = 1 — k 2z /w 2 . The k 3 and k± integrals are Lorentz invariant and can be 
evaluated in the centre-of-final- momentum frame, as in (|A14 ). Hence 

_ V r Vd 3 k 3 r Vd?k A 2 
° ~Tv~lJ J "(2^' 1 

\g,g u \ 2 cl f d% kje(kj) {kl _ k2)l/ M* (Bg) 



16vr 2 M 2 J \k] -m 2 +im^| 
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6 is the Heaviside step function. 
As Y^/m^ -> 



\kj - ml + im^T^ 2 T 



8(k 2 - ml) + 0{l) . (B9) 



If tuf — k] z < mi the first term of ( |B9|) does not contribute to (|Blj| ) and cr ~ \g u g u \ 2 / M 2 . If 

w 2 — kj z <C m 2 . then a ~ \g u g u \ 2 (wj — k 2 Iz ) 2 j '(M 2 mi). Otherwise, to leading order, the 0(1) 
terms can dropped and 

° = |^J 2 ^i^e-^ + ^)^ 2 / (2A ;7r A; 2r /M 2 ) (BIO) 

with kj T = wf n — kf nz — m 2 . The angular integral is evaluated using J 27r d6e xcos9 = 2irl (x). 
The plasmajparticles will have energies of order T, so if w\ is less than the GUT scale then 
kiT, &2T < M and the exponential terms can then be dropped. 

Scatterings involving charm and top quarks will give similar contributions to ( |B10| ). 
Summing these and substituting (|A16|) gives 



c 2 



a 



~± (©« - kL - ml) + 0(\g\ 2 ) + 0(T Wl /M 2 )) . (Bll) 



The above calculation only involves the neutral component of the intermediate Higgs 
field. Similar scatterings involving the other components also occur, such as u^ tr +e — > d+u c . 
In SO (10) there are a total of 5 such scatterings. They all have similar cross sections to 

(EH- 



It is also possible for plasma particles to interact with the neutrino currents by exchanging 
a virtual Higgs particle, as in figure |3|b. Conservation of momentum implies that the Higgs 
particle is always space-like, and so the scattering is never resonant. Thus the dominant 
contributions to the total cross section come from ( |B11| ) and similar interactions. 

The rate at which the currents interact with the plasma is given by r int = (crt> rc i), 
where is the equilibrium number density of each species of plasma particle and (cxf re i) 
is the thermally averaged cross section. n p i as ~ T 3 and a ~ \g u \ 2 /M 2 , so 

r int ~ ^£t 3 . (B12) 
M 2 1 ; 



APPENDIX C: DECAY OF A MASSIVE BOUND STATE ON A STRING 

Conservation of energy-momentum prevents massless bound fermions decaying. The 
same is not true of massive bound neutrinos. These can decay into left-handed neutrinos 
and light Higgs particles v^ tr — > v + $ (see figure ||c). Since they are Majorana fermions, 
they can also decay into the corresponding antiparticles z/ c tr — > v + $. Both amplitudes 
involve similar calculations. We will just consider the lowest mass bound state, with angular 
dependence n — — 1. Typically the heavy neutrino will have a vacuum mass of order the 



22 



GUT scale. The mass of a massive bound state will then be some fraction of the GUT 



scale. Using the expressions (|3.4j) , (|A3|) , ( JA8|) , the amplitudes for the above decays in the 
rest frame of the bound state are 

A = k c tr (k 1 )|5|p(k 2 )$(k 3 )) 



. J d'x C/J kl (^ia'V^^^e^^-^ 

J Vo™3^>' 



= %gl J d'xUl^i^U^^^e^-^-^ 



= igl I Sx ^ e^^)^ - * a - fe 8 ) Jl - ■ (C2) 

As with the massless current scattering (|B3| ), we will use a Bessel function expansion ( A17 ) 
of e i ( k2+k3 )' x(2) , and approximations of the radial integrals. 



J (A;r)x2(Msr)rdr « ^e" fc ^ j , y Ji(kr) X i,3( M s r ) rdr ~ ~ 2 ^, e~ fe /i2M } , (C3) 
with C2 ~ —0.7 and ci ~ c 3 ~ 0.3. Thus 



Ax - -^^ (2?r)3 \.UJ^ Cx hl±^3L + ./l+Scae^ 
MVy/8nw 3 L \ V w 2 M B V ^2 



x5 (M)( fcl _A; 2 -A;3)e- fe '/( 2M2 ) , (C4) 



^ ~ - ^^( 2?r ) 3 /1 _ S^ ^d^ ^M^ _ fc2 _ fc 3 )e-^ 2 /(2M 2 ) ( C5) 
MUv^^LV »2 M B 1123; , 

where {k s ) x ^ y = (k 2 + k 3 ) XjV . As with the massless scattering in appendix [B], transverse 
momentum is not conserved and there is a M _1 suppression due to the bound state being 
confined to the string core. The 2-component vector k s gives the momentum contributed 
by the string. The two amplitudes correspond to different decay products so the total decay 
rate is obtained by squaring and then adding them, and finally summing over the outgoing 
state momenta. Before doing this, it is useful to introduce 2 delta functions to separate out 
the k s dependence. 

1 fV 2 d 3 k 2 d 3 k 3 



T -tJ (2*)° U^l + \ A 

= > |2 [ d%^^6\h + k s -k 2 -k 
4M 2 (2vr) 3 J w 2 W3 

( k 2 k ■ k 



1/r ^ + c> 2 + 2c lC2 ^ + c> 2 )e-^ 2 . (C6) 
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The d 3 k 2 and d 3 k 3 integrals can be evaluated using 



d 3 ko d 3 k 



^±^± k »5\k - k 2 - k 3 ) = k»7r I 1 - -J ] Q(k 2 - ml) . 
w 2 w 3 \ k 2 I v 



2\ 2 



(C7) 



This can most easily be shown in the centre of outgoing momentum frame, and then Lorentz 
transformed to a general frame. We expand the remaining d 2 k s integral as a power series in 
m^/M B . 



\g„\ 2 M B fV M B~ m : 



16ttM 2 Jo 

12 



kgdkg 



+ (c? + c| + 2cic 2 



Mi 



-k 2 JM 2 



M 



B 



32tt 



{c|(l - e "*) - (c? + 4 + 2c 1 c 2 )[A- 1 + e"*(l - A" 1 )]} + O 



(C8) 



where A = (M B /M) 2 ~ 1. 

The massive bound state lifetime is then 



r ~ (\g u \ 2 M B 



(C9) 



which is small. 



APPENDIX D: ZERO MODE SCATTERING ON NEARBY STRINGS 

In addition to plasma interactions, it is also possible for massless fermion charge- carriers 
to decay by interacting with currents on other strings. This will be most significant when 
the separation of the strings is small. 

We will consider two infinite straight strings with opposite windings, running parallel to 
each other. We define 2R to be the displacement of the second string, with R orthogonal to 
the z-axis. This calculation will also be of relevance to string loops, R is then the radius of the 
corresponding loop. For a typical vorton (a loop stabilised by a current) R ~ 10 — lOOM^ 1 . 




Since the second string has opposite winding to those considered in section JTJ, the 
solution (|3.5|) is not valid. Instead 

"ok K ~7=t I i ] exp (-/£/ + (Dl) 

should be used, with p 2R a radial coordinate centred on 2R rather than the origin. These 
states have k z < 0, and opposite chirality to (|3.5| ). 

The cross section for direct interactions between currents on different strings will be pro- 
portional to the square of the overlap of their wavefunctions, which is of order exp(— AMR). 
Interactions involving the exchange of a virtual particle do not suffer from this suppression, 
and provide the dominant decay channel. We will consider the exchange of a virtual Higgs 
particle v c str + V^ tr — > v + i>, as shown in figure ^]d (the arrows indicate the direction of the 
bound state momenta along the strings). 
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A = (v:Mv:M\s\v(k 3 )v(k 4 )) 

= 2\g u \ 2 Y: J d A x J d'y (^(x^E^CxJe^^^CxOi^V^C^e^-^ 

k 

- ^ 1 (x)^V i r 4 (x')e i (---^S k2 (x')^ 2 f/ k3 (x)e^--^)*') G (x - x>) . (D2) 



Substituting (|3iS|), (|A3|) , (|A8| ) and (|13 1| ) , with [/ ok = f/ ok , and evaluating the integrals gives 



47r|^| 2 cg 




L\/M 2 





1 + — \/l - — - fe,, ^2 - ^ 3 )e 2jk3 - R ^ <^ 2) (^i + h - k 3 - h) (D3) 

w 3 V w 4 



where 



, f d 2 k f (k + k 3 )l (k-k 4 )l \ . . 

l(k z , w) = / — 9 -r exp - l \J )T - K - ^J^- - 2ik ■ R . D4 

y ' J m 2 + k 2 z - w 2 + k 2 , \ 2M 2 2M 2 / 

Applying the convolution theorem gives 

l(k z , w) = 2M 2 J d 2 XK (2^m 2 +k 2 -w 2 \R - X|) 

xe -X 2 M 2 e i(k 3 -k 4 )-X e -(k4+k 3 )|/4Al 2 _ / D5 x 

The main contribution from the integrand occurs around X = 0. Since the modified Bessel 
function K varies slowly in this region, ( P5|) can be approximated by putting X = in the 
argument of the Bessel function. This gives 



l{k z , w) w 2ne- k lTl 2M2 e - k lT/ 2M2 K (2 Jm 2 + k 2 - w 2 R) . (D6) 



The exponential factors come from non-conservation of transverse momentum at each string 
and the Bessel function is a result of the limited range of the virtual Higgs particle. 

Substituting (|D6l ) into (|D3|), squaring and multiplying by L/Tv re i gives the cross section. 
u re i = 2 and so, 



2| ^i? / dw 3 dn 3 dw A d^ z \h + k 2 -k 3 - h)w 3 w 4 e- k2 *T/ M2 e - k iT/M 2 
(2vr) 4 M 4 V 

{(w 3 - fc 3 *)(w 4 + A; 4z )Xo + (k lz - k 3z ) 2 - { Wl - w 3 ) 2 R) + (k 3 < — ► k 4 ) 



-2k 3T k 4T Re ( e ^-2k3-R) e ^ 4+ 2k 4 .R)) Kq + (k lz - k 3z ) 2 - (w, - w 3 ) 2 R) 

K {2^m 2 + (k lz - k 4z ) 2 - ( Wl - w 4 ) 2 R)} . (D7) 

The cross section gives the transition rate for unit incident flux. The incident flux in 
this case is L/v Te \, instead of V/v re \ as in ( [Blip , because the incoming particles are confined 
to the one- dimensional strings. This gives a dimensionless cross section instead of the usual 
(length) ~ 2 . 
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The Bessel function means that the first term of ( |D7| ) is exponentially suppressed except 
when (ki z — k 3z ) 2 — (wi—w 3 ) 2 is small. We will use the approximation K (2z) = — log(z) 
for z < 1 and K (2z) = elsewhere. In this region the transverse outgoing momenta are 
small, and the gaussian factors can be dropped to first order. If m^R <C 1, the Higgs mass 
can also be neglected. 

The integral of the second term of (|D7|) is equal to the integral of the first. The third 
is of order k\ T k\ T ^ which is small in the region where the Kq factors are not exponentially 
small. Thus the third term of ( |D7| ) can be neglected. Using the above approximations gives 



where u>coM is the incoming centre of mass energy. 

If m^R is large the above approximations do not apply. This occurs if $ is a GUT boson, 
or the electroweak Higgs field at high temperatures. In this case an asymptotic expansion 



Taken at face value, ( P8|) gives a short lifetime for charge carriers on a vorton. However, 
there are complications in directly applying this result derived for straight strings to curved 
vortons. The above calculation does not take into account conservation of this angular mo- 
mentum, as the system considered has no rotational symmetry. For a perfectly circular loop 
we would expect angular momentum to be conserved. Combined with energy conservation, 
this would prevent massless modes on the string scattering into massive modes. However, 
the fermionic spectrum has not been calculated for a circular loop and there is no reason to 
expect the zero modes to remain massless. Angular momentum conservation should only be 
considered in a consistent rotationally invariant calculation. The calculation above works 
consistently with straight strings, thus while ( P8|) may not be directly applicable to vorton 
decay, it is relevant for interactions on non-circular loops and scattering of currents on a 
string network. 




(D8) 



of the K factors can be used to estimate 




9 A m 4> c -4m A B 

M 4 i? 3 



(D9) 
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